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Motivated by cosmological considerations, and by the atmospheric and solar neutrino flux deficits, 
we consider the scenario in which the three Standard Model neutrinos are endowed with a nearly 
equal Majorana mass in the eV range. Combined constraints coming from (1) direct search for 
electron neutrino mass, (2) absence of neutrinoless double beta decay, and (3) unsuppressed electron 
neutrino flux in present reactor experiments, imply a quite well specified pattern of neutrino masses 
and mixing angles. We discuss the experimental tests of the model and comment on the features of 
this scenario. 



I. NEARLY DEGENERATE MAJORANA NEUTRINOS 

A. Phenomenological character of the approach 

The possibility that the neutrinos and u^- produced in weak interactions are coherent combinations of almost 

degenerate mass eigenstates has attracted great attention in the particle physics community p|-p^ . This is motivated 
by the fact that the mass splitting required to explain the atmospheric and the solar neutrino flux deficits in terms of 
neutrino oscillations arc much smaller than the mass required for cosmological purposes, as will be discussed below 
in this first section. 

Most of the works focused on the ansatze for the neutrino mass matrix: = mg x 1 + small corrections, and on 
the study of the corresponding theoretical motivations and phenomenological implications. Inevitably, this ansatze 
implies a signal for neutrinoless double beta decay {Qvj3(3) experiments if toq = few eV's, of interest for cosmological 
considerations. Allowing a departure from this mass pattern, the authors of ^ outlined a scenario in which the 
neutrinos have masses in the sub-eV range and could produce a signal in next generation searches of Qv(3l3. It was 
later proven that it is possible |jll| to incorporate a massive neutrino in the eV range in the context of models with 
abelian family symmetry, in a way marginally consistent with the 0vj3(3 bounds. 

In the present work we do not commit a priori ourselves to any given pattern for the neutrino mass matrix, and 
concentrate on finding a picture of the neutrino masses and mixing angles that is consistent with the phenomenological 
observations, and with the hypothesis of massive neutrinos in the cosmologically interesting range. We will show that 
these hypotheses are quite restrictive, implying that the corresponding scenario is likely to be confirmed or disproved 
by forthcoming experiments. 

In its phenomenological character, our approach is very similar to the one in [ p2[ . We emphasize in the present 
study the importance of the direct search for neutrino masses, as crucial test of the scenario under discussion. Another 
difference with is in the fact that the authors of that work concentrate on the MSW solution of the solar neutrino 
problem [p"3|-|l5|, whereas we focus on the vacuum oscillation solution [p^ ; in this respect, the two analysis complement 
each other. 



B. Definition of the scenario 



The lagrangian of neutrino masses in the basis where the charged lepton masses are diagonal is (bispinorial notation) : 

i:mass = - Mil. Ni,{x)Ni,,{x)+\l.C.] (1) 

£ and £'—e,fi,T 

the Lorcntz indices are contracted according to N^Na, A = 1,2. The symmetric matrix A4 can be decomposed as: 

Mw = ^i'^ (2) 

1=1,2,3 

(in compact notation: M. = U*mU^), where U is the unitary matrix of mixing angles and phases. The three 
nonnegative and nondecreasing real numbers m^. are the Majorana neutrino masses. Introducing the bispinors Ni{x): 
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1=1,2,3 



(3) 



and the Majorana four-spinor in the chiral representation; Viaix) = (A^iyi(x), Nf-{x)) we can finally rewrite the mass 
lagrangian in the standard form: 

-Cmass = ^ in^.V-{x) iyi{x). (4) 

1=1,2,3 

The mass difference is denoted as: 

Ato| = ml^ - ml_ . (5) 

It enters for instance in the expression P«/ = | C/^* [/^/j exp(— of the probability of oscillation in vacuum 
from a flavour £ to £' . In fact we can replace: exp{—iEit)/ cxp{—iEjt) by exp{—iAmfjt/{2p)), under the assumption 
(which holds in the cases of interest) that the neutrino momentum p is much larger than the neutrino masses. 
The scenario under discussion can now be defined by letting in first approximation: 

TO I,. = mo z = l,2,3. (6) 

The value of too is fixed below, with the help of cosmological considerations. Before passing to that, we further specify 
two sub-cases of this scenario, by making the assumption that the squared mass differences are in ranges of interest 
to the solar and the atmospheric neutrino problems. Case A is defined by: 

Ato^i = Ato| « 10-i°eV2, Aml^ = Aml,^ « IQ-^eV^; (7) 

in case B the two mass differences are inverted (compare with fig. 1 in last section). Since, apart from this fact (and 
neglecting matter effects on the propagation of neutrinos), the two cases are phenomenologically equivalent, we will 
concentrate the discussion on case A. 



C. Neutrino masses of cosmological interest 

There are various indications of the existence of dark matter (DM — see for instance [0), that may constitute the 
main component of the Universe. Once the amount of dark matter, ordinary matter and, possible other components 
are specified, one can attempt a comparison of the observed structure with the one predicted in terms of the spectrum 
of the primordial fluctuations. The gross features of the Universe can be explained by the cold DM model, that 
assumes (1) that the dominant component of DM is cold (non relativistic at the moment of decoupling), (2) a nearly 
scale invariant (Harrison-Zel'dovich) primordial spectrum, (3) and that the parameter is 1, which implies that 
the total energy density is equal to the critical density pc (the latter two hypothesis are suggested by inflationary 
scenarios). However COBE data have indicated the need of a certain departure from this simple cosmological model 
I p^ . The hot+cold DM model jl^,^ achieves an optimal agreement with the data if the hot DM constitutes the 
20 — 30% of the energy density, the cold one the 70 — 60%, and the baryonic matter approximatively 1 tenth of the 
total energy density^. 

Massive neutrinos are natural candidates for hot DM, and this is the assumption we are interested to consider. The 
number density of neutrinos can be computed in the Big-Bang model p^ . Once they become non-relativistic, 
their energy density reduces simply to pi, = TOj/.rij^. Assuming Q — I, and introducing the ratio of to the critical 
density pc (proportional to Hq) one arrives at the expression that relates the neutrino masses to the cosmological 
parameters: 

m^. — 6.6 eV x 



0.2 Pc 



0.6 



(8) 



'^In this model, the favoured value of the Hubble constant Hq = hx 100 km sec~^ Mpc~^ turns out to be around h = 0.5. In 
pl| the observational range h = 0.65 ~ 0.85 is suggested; a more recent analysis |Q gives h = 0.6 — 0.8 as the reference range. 
Even lower values, h — 0.5 — 0.7 have been advocated in [ESll. 
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The neutrino mass differences implied by solar and atmospheric neutrino deficits are small in comparison with 
the mass in eq. (||), and this facts suggest the hypothesis of nearly degenerate neutrinos. Using eqs. and (||) 
together with the informations from the study of the structure evolution, we are lead to the following reference value 
for mo = J2t "tiy,/3 : 

mo = TO(cosm.) — 2.2 eV. (9) 

This corresponds to h — 0.6 and — 0.2 p^, or, with approximation of 5%, to h = 0.5 and — 0.3 Pc. Progress 
in the determination of Hq and il will provide tests of the hot+cold DM scenario, and also, if the scenario will be 
confirmed, a better determination of the parameter mg. 

We briefly comment on the possibility that massive neutrinos may constitute the major component of the DM. It 
was suggested that this could happen if special seeds of primordial fluctuations like those originated by the cosmic 
strings ||25|| play a major role in the process of structure evolution, or, alternatively, if neutrinos participate in exotic 
long-range (nongravitational) interactions ]2^ ]. However, apart from any theoretical justification, if this alternative 
scenario is correct, and fl — 1, the reference value of the parameter mo would be larger than 10 eV, since in eq. 
increases almost 5 times. 

Of course, it is a major assumption to consider the absence of light sterile neutrinos, that could play various roles, 
and in particular could give the dominant contribution to the hot DM content of the Universe [p8|-p0[ . The hypothesis 
of (Standard Model) singlet neutrinos could be tested if they are mixed with the ordinary ones, and participate in the 
oscillations. Regarding this point, if neutrino masses are the explanation not only of the solar and of the atmospheric 
neutrino problems, but also of the LSND data on v^j^ — Vf. [ ^T| , ^ and on v^^Ve flux conversion psf , major modifications 
of the scenario, like the introduction of a sterile neutrino would be required. Therefore, in order to check the scenario 
under discussion (in which no light sterile neutrino is present and a negligible signal at LSND is expected) it will be 
crucial to see if the LSND findings will be confirmed by the forthcoming upgraded KARMEN and/or by other 
experiments. 

II. CONSTRAINTS FROM PHENOMENOLOGY 

In this section we discuss the bounds on the neutrino masses and mixing angles that can be deduced from the 
absence (1) of positive indications of electron neutrino mass in the tritium decay experiments; (2) of neutrinoless 
double beta decay signal; (3) of appreciable disappearance of electron neutrino flux in reactor experiments. The 
constraints are considered in order of relevance for the scenario with three almost degenerate neutrinos^. 

If one (or more) of the related experiments finds a signal, one (or more) parameters of the model will be determined. 
We stress that, in the model under discussion, (i) the direct detection of neutrino mass seems to be possible, and 
that {ii) condition (2) implies a strong restriction on the mixing angles, that will further strengthen if future Oi^/3/3 
experiments will continue to give a null result. 

A. Direct search for neutrino masses 

Let us denote by dr[my^]/ dEf. the differential electron spectrum in beta decay supposing that Vi coincides with the 
electron neutrino. In the general case, in which the electron neutrino is not a mass eigenstate, the beta spectrum 
dV/dEe is the weighted combination: 

^ = E - - ^e)' (10) 

i=l,2,3 

where the Heaveside function Q{x) enters the formula due to kinematical requirements, and AA/ is the difference 
between the mass of the decaying nucleus and that of the final nucleus. In our case, due to (||), the unitarity relation 
X^i Wei\^ = 1 implies a constraint straightforwardly on the common mass mo. 

Stringent limits can be derived by the experimental study of endpoint spectrum in tritium decay, after taking into 
account in some way the problematic indication of negative masses squared common to most recent experiments 



^The importance of the direct searches for neutrino mass has been not, in our opinion, sufficiently emphasized in the literature, 
even if the existence of this constraint was realized since the beginning (see for instance fig. 2 in reference Q). 
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(corresponding to the presence of a bump in the endpoint region). In the most recent experimental analysis 
still under progress, new parameters unrelated to the neutrino mass are introduced^ The upper limit obtained, 
ml — 1.5 ± 5.9(stat.) ± 3.6(syst.) eV^, implies: 

mo < m(3H) = 3.9 eV at 95% CL, (11) 



summing the errors in quadrature. The limit (11) is not very different from those obtained by other experiments; it 
could be further strengthened to 2.2 eV if the bump would be fully understood ||35[] . 

The limits on neutrino mass from SN1987A obtained by the measured spread in time (which is related to the 
original spread in energy and to the duration of the neutrino pulse) are significatively looser, to(SN1987A) ^ 15 eV. 
However, the observation of neutrinos from a very distant supernova could yield a competitive or even better bound 
in future. 

The following remarks should clearly illustrate that, in the scenario with nearly degenerate neutrinos, the direct 
search for neutrino mass is of cosmological relevance: 

1. The "preferred" value (^) is quite close to the limit (pi]): mo « 2 x m(^H). Therefore, the hot+cold DM model 
can be ruled out or confirmed by a slight increase in the precision of the measurements. We stress once more the 
importance of a clear theoretical description of the structure of the endpoint spectrum in tritium decay studies. 



2. For pi, = 0.2 pc the limit (|ll|) disfavours h > 0.8 in the hot+cold DM model. Increasingly precise measurements 
of the cosmological paramenters and numerical simulations will be important in narrowing the parameter space 
of the model. 

3. The scenarios in which pi, dominates the energy density of Universe and that have $7 = 1 require toq ^10 eV, 
and therefore can already be excluded in the present context. 

B. Neutrinoless double beta decay and large mixing 

The limit (|l^) (the preferred value (|^)) are approximatively 6 times (3 times) larger than the limit obtained by 
searches of neutrinoless double beta decay p7[|: 



< m(0z//?/3) = 0.68 eV at 95% CL. (12) 



We disregard in the following the possibility of large errors in the estimation of the hadronic matrix elements. The 



bound (12) implies a cancellation in the contributions to |Alee| with a precision of ■m{Qv (3 (3) / ttiq ~ 30 %. Let us first 
analyze the case in which the neutrinoless double beta decay is strongly suppressed, which implies C/gj ~ 0. This 
equation can be represented geometrically as a triangle in the complex plane. Therefore it can be satisfied only if the 
sum of the length of two sides is not smaller than the length of the third side. Taking into account that unitarity fixes 
to unity the perimeter of the triangle, we conclude that the values of the moduli for which cancellation is possible 
must satisfy: 

\>\U.i\\ 

1 1 (^^^ 

Assuming instead that a signal is observed, the upper bound on the elements of the mixing matrix resulting from 
geometrical considerations reads: 

|c/.,r<i±^^:^i^. (14) 

Even with a maximal value of the mass ratio: |Alee|/™o = 30 % (that would imply that a signal of Ovjdfi will be 
observed close to the present limit) the upper bound does not change much: |C/ejP < 0.65. 



^In last analysis these parameters have be related to instrumental or maybe new physics effects. 



4 



We want to make one point clear: we are not proposing any explanations of the cancellation between the elements of 
the mixing matrix squared, we are just remarking that, once the hypotheses underlined above are accepted, consistency 
with the data requires a cancellation; this is true even if a positive signal, related to a small mass A4ee should be 
found. We show below that quite a precise pattern of the mixing matrix can be figured out insisting on this scheme, 
and on the consistency with the experimental data. 

Eq. (^3|) implies that there is at most one small mixing matrix element Uej , and even in this case the other two 
elements correspond to maximal mixing angle. The presence of large mixings, that was generally noticed in previous 
studies is simply a consequence of the absence of signal in Oi/Pf] searches in these models. 

It was remarked that the scenarios with large lepton mixing tend to be disfavoured by the comparison of the 
observed SN1987A neutrino energy spectrum and the theoretical spectra [Q. Therefore, a clear understanding of the 
mechanism of neutrinos emission and propagation in the process of stellar collapse could provide us with the most 
direct way to exclude the neutrino scenario under discussion. We will not speculate further on this idea in the present 
work, and concentrate instead on the study of possible patterns of the mixing matrix. 

C. Information from the reactors 

The neutrino beams produced at nuclear reactors are well suited for the study of disappearance of electron neutrino 
flux. The neutrino oscillations are described by the conversion probability 



Pee = l-4^|C/e,|2|t/e,:|2sin2 

j>i 



4 p 



(15) 



which just depend on the mixing angles which are constrained by the neutrinoless double beta decay. In the case A, 
since that the parameter AtUq is too small to develop oscillating phases on terrestrial distances, our expression ( p^ ) 
reduces to: 



Pee = l-4|[/e3p(l- If/esHsin^ 



4 p 



(case A) (16) 



that is effectively a two-flavour neutrino task, the relevant element of the mixing matrix being the one related with 
the "isolated" massive neutrino (in case B, one arrive at the same expression with Ues replaced by Uei)- The neutrino 
flux is not suppressed at all in the two extreme cases: |?7e3p « or « 1, but the second solution must be discarded 
in view of the absence of OvP/S signal, that implies eq. (|l^). Applying the inequalities ( |l3| ) we conclude: 

|?7eip«l/2, |C/e2p«l/2, |C/e3p«0 (casc A) (17) 

(in case B the small element is |C/ei|)- The constraints from reactor experiments are not, of course, absolute. To be 
quantitative we notice that for a value of Arn^j^ = 0.01 eV^ that may explain well the atmospheric anomaly we get 
the bound [|||: 

|J7e3p < 0.2 at 90% CL (case A). (18) 

This bound however depends strongly on the value of Am^^^, and becomes already ineffective if instead of 0.01 we 
have 7.5 x 10^'^ eV^ : This should warn from a literal use of equation ( |l7|) (compare also with the discussion in |p^). 
The future search at Palo Verde and at CHOOZETj) will be able to strengthen the bound, or maybe suggest 
a non-zero value of the mixing angle. Assuming eq. (|lj)7thc bounds from neutrinoless double beta decay are easily 
solved: Ue2 ~ ±z x Uei, where i is the imaginary unity. 

III. FEATURES OF THE SCENARIO 

Let us summarize the major consequences of the constraints exploited in previous section by using the matrix \U\, 
whose elements are the numbers \Uej\- Using eq. (p^, assuming for simplicity that Ues is quite small — our stronger 
hypothesis — , and the unitarity relationfQ the approximate pattern is: 



*In strict analogy with the what happens to the Cabibbo-Kobayashi-Maskawa mixing matrix, the matrix \U\ can be described 
in the general case by four parameters (for instance one can use three mixing angles and one phase). 
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\u\ 



I s/2 V2 



sm 9 

cos COS Q 
\ V2 V2 



V2„ 




(case A) 



(19) 



(in case B, the last column exchanges with the first one). There is only one parameter in eq. ([l9|), Q £ [0, 7r/2]; we 
will see in the following discussion that the preferred value is « 7r/4. As discussed in sub-section C. above, the zero 
entry should be intended as an approximate statement only. In the rest of this work we will discuss the qualitative 
features of the scenario outlined, leaving the detailed comparison with the data for a future work. 

A. Atmospheric neutrinos 

On terrestrial distances L w _Re , and for neutrino momenta p « 1 GeV the oscillations related to the mass difference 
Am0 do not develop: cos(AmQL/(2p)) 1. Therefore, in case A, the expressions for the matrix V of the conversion 
probabilities (whose elements are where £ and I' — e,/i, r) is: 



V Ril-sin^ 26* X sin^ 



4p 




(20) 



If we suppose that the oscillations related to Am^^^ are effectively averaged, at least at distances L of the order of 
the Earth size, we can replace the oscillating factor sin^(Am^t„L/(4p)) by 1/2. A rough estimation of the muon-to- 
electron neutrino ratio r = $(i'^)/$(z^e) measured in the detectors can be obtained assuming that in the production 
region rg = $o(t'/i)/^o(^'e) ~ 2, that is adequate for neutrinos of energies around 1 GeV originated by cosmic rays 
entering the Earth atmosphere. Using ^{yi) — ^ P«' $o(^'£') we find for the double-ratio R : 



R 



£'—e.fi,r 

_ r 



1 



sin2(26l) 



(21) 



Kamiokande and Superkamiokande ||4^ measurements in the multi-GeV energy region suggest R ~ 0.55, and this 
implies that 9 has to be close to = 7r/4 to produce a large suppression factor. However, the data are much better 
fitted if a linear dependence of the spectrum from the cosine of the zenith angle (j>z is also introduced: 



R{(j)^) « 0.55 + 0.3 X cos 02 



(22) 



This favours not-so-large values of Am^j^, in such a way that the averaging effect is reduced for smaller neutrino 
travel distance (around 10 — 100 km, corresponding to 4>z in a neighbourhood of 0) and it is maximal for neutrinos 
produced in the antipodal regions (produced at L ~ 10000 km far apart from the detector, that is for (j)z ~ tt). So 
the averaged oscillations would take place at maximal distances only: R{n) « 0.3, which should be compared with 
the smallest possible theoretical value, R ~ 0.5. This is not excluded by the data, the errors in each of the 5x2 data 
points being around 0.1 (except for (j^z ~ 0, where the errors are 4 — 5 times larger). But while waiting for more 
precise data, one has to keep in mind that the model considered suggests somewhat smaller flux suppression than 
measured at present. 



B. Solar neutrinos 



We consider the vacuum oscillation mechanism as a solution of the solar neutrino problem, which requires Amg 
10^^*^ eV^. The survival probability of electron neutrinos in case A is: 



Pee = l- 2|C/e3p(l - |C/e3p) - 4 1 C/el P | [/e2 ^ siu^ 



Am|I/ 



4 p 



(23) 



where L represents the Earth-Sun distance. The second term gives the oscillating factor, whereas the first term can 
give an additional (constant in energy) suppression. In fact, if |f/eiP — |C^e2p we can rewrite: 



Pee-|(7e3|^ + (l-|[/e3n'cOs2 



Ato|L 
4 p 



(24) 
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If the bound ( Jl&j ) holds, the first term is neghgible; therefore the oscillating factor can be multiplied by a suppression 
factor that can be as small as 0.64. This can somewhat diminish the spectrum distortion and the seasonal variations 
that are correlated characteristics of the vacuum-oscillation mechanism [Q. If, beside that, we suppose that the 
oscillations are effectively averaged (requiring that the splitting Amg is large in comparison with 10~^° eV^), the 
probability of suppression is close to the maximal averaged value {Pee)max = 1/3. This may be of relevance for 
neutrinos emitted by a supernova, or maybe for solar neutrinos themselves p6| . 



IV. SUMMARY AND DISCUSSION 



We summarize our findings in fig. 1, which visualizes the flavour content \Uii\'^ of the massive neutrinos Vi in the 
way used in to present a panorama of scenarios of neutrino masses and mixing. 



ml. — ml (case A) 



ml. — ml (case B) 





1^3 



^1 



1^3 



FIG. 1. Composition of massive neutrinos: white for electron-, hatched for muon- and black for tau-neutrino flavour. For 
a proper interpretation of the vertical scale, recall that the mass splittings, eqs. (^) are much smaller than the common mass 
scale mo, eq. (|9|). 

Let us recall the main features of the scenario proposed, that assumes that the solar and the atmospheric anomalies 
in the fluxes of neutrinos are due to oscillations: 

1. Cosmological indications can be reconciled with the the solar and atmospheric neutrino deficits if the three 
Standard Model neutrinos have in first approximation the same mass toq = m(cosm.). 

2. We emphasized that the common mass mo needed in hot+cold DM scenarios can be probed in direct searches 
for neutrino masses, and that this mass is only two times smaller than the masses that can be currently excluded 
(eqs. (|) and (|ll|)). The case in which the energy of the Universe is dominated by three massive neutrinos can 
already be excluded. 

3. The absence of a signal in neutrinoless double beta decay searches requires a cancellation. Even if this is not a 
priori a theoretical attractive perspective, it implies strong restrictions on certain combinations of the elements 
of the mixing matrix (compare with eq. (p^). 

4. The absence of signal in reactor experiments suggests that the structure of the mixing matrix, eq. (^9|), may be 
quite simple (see also eq. ( |l8[ ) and related discussion). Regarding this point, future electron neutrino disappear- 
ance experiments will be important to better fix the parameters of the model. Even better, they could permit 
us to distinguish the present variant of the scenario with nearly degenerate neutrinos with the variant in which 
the solar neutrino deficit is explained invoking the MSW effect; in fact, in this latter case large values of the 
mixing element Ues are preferred (compare with eq. ( p^ ) and discussion thereafter) . 

5. The pattern of the mixing matrix allows us to take into account the atmospheric neutrino anomaly. With 
accumulation of data at SuperKamiokande, the zenith angle dependence of the neutrino spectrum will become 
a severe test of the scenario proposed. 

6. Regarding the observed suppression of the solar neutrinos flux, we remarked that this model can incorporate 
the vacuum oscillation solution, possibly with reduced distortion of the energy spectrum. 
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In conclusion, we express our opinion that, even if the model proposed is not at present of particular theoretical 
appeal, it is of some interest since (a) it incorporates cosmological and phenomenological indications; (b) it is enough 
restrictive to be confirmed or disproved in the future (few) years. 
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